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Abstract 



These lectures are designed to provide a general introduction to the Einstein-Vlasov 
system and to the global Cauchy problem for these equations. To start with some general 
facts are collected and a local existence theorem for the Cauchy problem stated. Next the 
case of spherically symmetric asymptotically flat solutions is examined in detail. The ap- 
proach taken, using maximal-isotropic coordinates, is new. It is shown that if a singularity 
occurs in the time evolution of spherically symmetric initial data, the flrst singularity (as 
measured by a maximal time coordinate) occurs at the centre. Then it is shown that for 
small initial data the solution exists globally in time and is geodesically complete. Finally, 
the proof of the general local existence theorem is sketched. This is intended to be an 
informal introduction to some of the ideas which are important in proving such theorems 
rather than a formal proof. 
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1. Introduction 

The Vlasov equation arises in kinetic theory. It gives a statistical description of a 
collection of particles. It is distinguished from other equations of kinetic theory by the fact 
that there is no direct interaction between particles. In particular, no collisions are included 
in the model. Each particle is acted on only by fields which are generated collectively by all 
particles together. The fields which are taken into account depend on the physical situation 
being modelled. In plasma physics, where this equation is very important, the interaction 
is electromagnetic and the fields are described either by the Maxwell equations or, in a 
quasi-static approximation, by the Poisson equation [26]. In gravitational physics, which 
is the subject of the following, the fields are described by the Einstein equations or, in 
the Newtonian approximation, by the Poisson equation. (There is a sign difference in the 
Poisson equation in comparison with the electromagnetic case due to the replacement of a 
repulsive by an attractive force.) The best known applications of the Vlasov equation to 
self-gravitating systems are to stellar dynamics [3]. It can also be applied to cosmology. In 
the first case the systems considered are galaxies or parts of galaxies where there is not too 
much dust or gas which would require a hydrodynamical treatment. Possible applications 
are to globular clusters, elliptical galaxies and the central bulge of spiral galaxies. The 
'particles' in all these cases are stars. In the cosmological case they are galaxies or even 
clusters of galaxies. The fact that they are modelled as particles reflects the fact that their 
internal structure is believed to be irrelevant for the dynamics of the system as a whole. 

These lectures are concerned not with the above physical applications but with some 
basic mathematical aspects of the Einstein- Vlasov system. First the definition and general 
mathematical properties of this system of partial differential equations arc discussed and 
then the Cauchy problem for this system is formulated. The central theme in what follows is 
the global Cauchy problem, where 'global' means global in time. Up to now global results 
have only been obtained in very special cases and one of these, the case of spherically 
symmetric asymptotically ffat solutions, is discussed in detail. The global existence results 
presented depend on having a suitable local existence result. At first the necessary local 
existence theorem is stated without proof and used in obtaining global theorems. The proof 
of the local existence theorem is discussed in some detail afterwards. Further information 
on kinetic theory in general relativity may be found in [10]. 

Let (M, (/o;^) be a spacetime, i.e. M is a four- dimensional manifold and g^p is a metric 
of Lorentz signature (— , -|-, -|-, -|-). Note that Qap denotes a geometric object here and not 
the components of the geometric object in a particular coordinate system. In other words 
the indices are abstract indices. (See [28], section 2.4 for a discussion of this notation.) It 
is always assumed that the metric is time-orientable, i.e. that the two halves of the light 
cone at each point of M can be labelled past and future in a way which varies continuously 
from point to point. With this global direction of time, it is possible to distinguish between 
future-pointing and past-pointing timelike vectors. The worldline of a particle of non-zero 
rest mass m is a timelike curve in spacetime. The unit future-pointing tangent vector to 
this curve is the 4- velocity f " of the particle. Its 4- momentum is given by mv'^. There 
are different variants of the Vlasov equation depending on the assumptions made. Here it 
is assumed that all particles have the same mass m but it would also be possible to allow 
a continuous range of masses. When all the masses are equal, units can be chosen so that 
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m = 1 and no distinction need be made between 4- velocity and 4-momentum. There is also 
the possibility of considering massless particles, whose wordlines arc mill curves. In the 
case m = 1 the possible values of the four-momentum are precisely all future-pointing unit 
timelike vectors. These form a hypersurface P in the tangent bundle TM called the mass 
shell. The distribution function /, which represents the density of particles with given 
spacetime position and four-momentum, is a non-negative real-valued function on P. A 
basic postulate in general relativity is that a free particle travels along a geodesic. Consider 
a future-directed timelike geodesic parametrized by proper time. Then its tangent vector 
at any time is future-pointing unit timelike. Thus this geodesic has a natural lift to a curve 
on P, by taking its position and tangent vector together. This defines a flow on P. Denote 
the vector field which generates this fiow by X . (This vector field is what is sometimes 
called the geodesic spray in the mathematics literature.) The condition that / represents 
the distribution of a collection of particles moving freely in the given spacetime is that it 
should be constant along the fiow, i.e. that Xf — 0. This equation is the Vlasov equation, 
sometimes also known as the Liouville or coUisionless Boltzmann equation. 

To get an explicit expression for the Vlasov equation, it is necessary to introduce local 
coordinates on the mass shell. In the following local coordinates on spacetime are always 
chosen such that the hypersurfaces x'^=const. are spacelike. (Greek and Roman indices 
take the values 0,1,2,3 and 1,2,3 respectively.) Intuitively this means that x^, which 
may also be denoted by t, is a time coordinate and that the x°' are spatial coordinates. A 
timelike vector is futiuc-pointing if and only if its zero component in a coordinate system 
of this type is positive. It is not assumed that the vector d/dx^ is timelike. One way of 
defining local coordinates on P is to take the spacetime coordinates x" together with the 
spatial components of the four-momentum in these coordinates. Then the explicit form 
of the Vlasov equation is: 

df/dt + ipyp'')df/dx'' - V'^ypWf/dp'' = (1.1) 

where are the Christoffel symbols associated to the metric Qafs- Here it is understood 
that is to be expressed in terms of and the metric using the relation QafSP'^P^ — 
An alternative way of coordinatizing the mass shell which is often useful is to use the 
components of the four-momentum in an orthormal frame, which has a priori nothing to 
do with the frame defined by the coordinates. It should be chosen so that the first vector 
is future-pointing timelike. The explicit form of the Vlasov equation in these coordinates 
is similar to (1.1), with the Christoffel symbols replaced by the Ricci rotation coefficients 
of the frame. Explicitly, it is given by: 

df/dt + {vyv^)e^df/dx'' - ^i^v^v^df /dv' = (1.2) 

The convention is used that frame indices are denoted by letters from the middle of the 
alphabet, while coordinate indices are taken from the beginning of the alphabet. The 
components of the frame vectors are denoted by and p'^ and are related by the 

equation p"' = e'^v^^. Since the frame is orthonormal v° = ^/l+~5^^^ViVj , where 5'^^ is the 
Kronecker delta. 
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The Vlasov equation can be coupled to the Einstein equations as follows, giving rise 
to the Einstein- Vlasov system. The unknowns are a 4-manifold M, a (time orientable) 
Lorentz metric g^is on M and a non-negative real-valued function / on the mass shell 
defined by Qce/s- The field equations consist of the Vlasov equation defined by the metric 
gai3 for / and the Einstein equation Gai3 = SnTap. (Units are chosen here so that the 
speed of light and the gravitational constant both have the numerical value unity.) To 
obtain a complete system of equations it remains to define T^p in terms of / and Qap- It 
is defined as an integral over the part of the mass shell over a given spacetime point with 
respect to a measure which will now be defined. The metric at a given point of spacetime 
defines in a tautological way a metric on the tangent space at that point. The part of the 
mass shell over that point is a submanifold of the tangent space and as such has an induced 
metric, which is Ricmannian. The associated measure is the one which we are seeking. 
It is evidently invariant under Lorentz transformations of the tangent space, a fact which 
may be used to simplify computations in concrete situations. In the coordinates (a;",p") 
on P the explicit form of the energy-momentum tensor is: 



Taf3 = - J fpaPp\g\^^^/podp^dp^dp^ 



(1.3) 



A simple computation in normal coordinates based at a given point shows that T^p defined 
by (1.3) is divergence- free, independently of the Einstein equations being satisfied. This is 
of course a necessary compatibility condition in order for the Einstein- Vlasov system to be 
a reasonable set of equations. Another important quantity is the particle current density, 
defined by: 

^" = -J fp''\9\^^Vpodp'dp''dp^ (1.4) 

A computation in normal coordinates shows that V^A^" = 0. This equation is an ex- 
pression of the conservation of the number of particles. There are some inequalities which 
follow immediately from the definitions (1.3) and (1.4). Firstly N^tV^ < for any future- 
pointing timelike or null vector V^, with equality only if / = at the given point. Hence 
unless there are no particles at some point, the vector A" is future-pointing timelike. 
Next, if V"- and are any two future-pointing timelike vectors then TapV^W^ > 0. 
This is the dominant energy condition ([12], p. 91). Finally, if X"' is a spacelike vector 
then Ta/jA^X^ > 0. This is the non-negative pressures condition. This condition, the 
dominant energy condition and the Einstein equations together imply that the Ricci tensor 
satisfies the inequality RapV'^V^ > for any timelike vector V". The last inequality is 
called the strong energy condition. These inequalities constitute one of the reasons which 
mean that the Vlasov equation defines a well-behaved matter model in general relativity. 
However this is not the only reason. A perfect fluid with a reasonable equation of state or 
matter described by the Boltzmann equation also have energy-momentum tensors which 
satisfy these inequalities. 

The Vlasov equation in a fixed spacetime is a linear hyperbolic equation for a scalar 
function and hence solving it is equivalent to solving the equations for its characteristics. 
In coordinate components these are: 

dX'^/ds = P" 

(1-5) 

dP^/ds = -F^ P^P^ 
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Let X°'{s,x°^,p°'), P^{s,x°',p°') be the unique solution of (1.5) with initial conditions 
= and P°'{t,x°',p"') = p°'. Then the solution of the Vlasov equation 
can be written as: 



where /o is the restriction of / to the hypersurface t — 0. This function /o serves as initial 
datum for the Vlasov equation. It follows immediately from this that if /o is bounded 
by some constant C, the same is true of /. This obvious but important property of the 
solutions of the Vlasov equation is used frequently without comment in the study of this 
equation. 

The above calculations involving T^/s and N°' were only formal. In order that they 
have a precise meaning it is necessary to impose some fall-off in the momentum variables 
on / so that the integrals occurring exist. The simplest condition to impose is that / has 
compact support for each fixed t. This property holds if the initial datum /q has compact 
support and if each hypersurface t = to is a Cauchy hypersurface. For by the definition 
of a Cauchy hypersurface, each timelike curve which starts at t = hits the hypersurface 
t = to at a unique point. Hence the geodesic flow defines a continuous mapping from the 
part of the mass shell over the initial hypersurface t = to the part over the hypersurface 
t = to- The support of f{to), the restriction of / to the hypersurface t = to is the image 
of the support of /q under this continuous mapping and so is compact. Let P{t) be the 
supremum of the values of attained on the support of f{t). It turns out that in many 
cases controlling the solution of the Vlasov equation coupled to some field equation in the 
case of compactly supported initial data for the distribution function can be reduced to 
obtaining a bound for P{t). An example of this is given below. 

The data in the Cauchy problem for the Einstein equations coupled to any matter 
source consist of the induced metric gab on the initial hypersurface, the second fundamental 
form kab of this hypersurface and some matter data. In fact these objects should be 
thought of as objects on an abstract 3-dimensional manifold S. Thus the data consist of 
a Riemannian metric gab, a symmetric tensor kab and appropriate matter data, all defined 
intrinsically on S. The nature of the initial data for the matter will now be examined 
in the case of the Einstein- Vlasov system. It is not quite obvious what to do, since the 
distribution function / is defined on the mass shell and so the obvious choice of initial data, 
namely the restriction of / to the initial hypersurface, is not appropriate. For it is defined 
on the part of the mass shell over the initial hypersurface and this is not intrinsic to S. 
This difficulty can be overcome as follows. Let be the mapping which sends a point of the 
mass shell over the initial hypersurface to its orthogonal projection onto the tangent space 
to the initial hypersurface. The map </> is a diffeomorphism. The abstract initial datum /o 
for / is taken to be a function on the tangent bundle of S. The initial condition imposed is 
that the restriction of / to the part of the mass shell over the initial hypersurface should be 
equal to /o composed with 0. An initial data set for the Einstein equations must satisfy the 
constraints and in order that the definition of an abstract initial data set for the Einstein 
equations be adequate it is necessary that the constraints be expressible purely in terms 
of the abstract initial data. The constraint equations are: 



(1.6) 



R - kabk""^ + (trA;)^ = 167rp 
Va/e^ - Vb{trk) = STTjb 



(1.7) 
(1.8) 



5 



Here R denotes the scalar curvature of the metric gab- If denotes the future-pointing unit 
normal vector to the initial hypersurface and /i"^ = g'^^+n'^n^ is the orthogonal projection 
onto the tangent space to the initial hypersurface then p = Tq,^?i"?i^ and = —h^^^Tp^n'^. 
The vector satisfies j'^Ua — and so can be naturally identified with a vector intrinsic 
to the initial hypersurface, denoted here by What needs to be done is to express p and 
ja in terms of the intrinsic initial data. They are given by the following expressions: 

P=j /o(p'^)p>a/(l +p''Pay/'{^'^9y/WdpV (1.9) 
ja = J fo{p'')Pa{^'^9)'/^dpVdp' (1.10) 

If a three-dimensional manifold on which an initial data set for the Einstein- Vlasov system 
is defined is mapped into a spacetime by an embedding iIj then the embedding is said to 
induce the given initial data on S if the induced metric and second fundamental form of 
il^{S) coincide with the results of transporting gab and kab with and the relation / = /oo0 
holds, as above. A form of the local existence and uniqueness theorem can now be stated. 
This will only be done for the case of smooth (i.e. infinitely differentiable) initial data 
although versions of the theorem exist for data of finite differentiability. 

Theorem 1.1 Let 5' be a 3-dimensional manifold, gab a smooth Riemannian metric on S, 
kab a smooth symmetric tensor on 5' and fo a smooth non-negative function of compact 
support on the tangent bundle TS of S. Suppose further that these objects satisfy the 
constraint equations (1.7)-(1.8). Then there exists a smooth spacetime (M, ga/s), a smooth 
distribution function / on the mass shell of this spacetime and a smooth embedding ip 
of S into M which induces the given initial data on S such that g^p and / satisfy the 
Einstein- Vlasov system and i^{S) is a Cauchy hypersurface. Moreover, given any other 
spacetime (M', g'^f^), distribution function /' and embedding iIj' satisfying these conditions, 
there exists a diffeomorphism x from an open neighbourhood of il^{S) in M to an open 
neighbourhood of ip'{S) in M' which satisfies x^fp = '4'' carries g^p and / to g'^f^ and 
/' respectively. 

The formal statement of this theorem is rather complicated, but its essential meaning 
is as follows. Given an initial data set (satisfying the constraints) there exists a corre- 
sponding solution of the Einstein- Vlasov system and this solution is locally unique up to 
diffeomorphism. There also exists a global uniqueness statement which uses the notion 
of the maximal Cauchy development of an initial data set, but this is not required in the 
following. The first proof of a theorem of this kind for the Einstein- Vlasov system is due 
to Choquet-Bruhat [7]. 

In the following we are mainly concerned with asymptotically flat spacetimes. These 
are the spacetimes which are appropriate for describing isolated systems in general relativ- 
ity. It is assumed that these spacetimes admit a Cauchy hypersurface with topology R^, 
although more general cases could also be considered. The smooth data set {gab, kab, fo) 
on is said to be asymptotically flat if there exist global coordinates x"' such that as |a;| 
tends to inflnity the components gab in these coordinates tend to Sab, the components kab 
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tend to zero, /o has compact support and certain norms are finite. These are weighted 
Sobolev norms. If tt is a smooth function on define 



m\H? 



"t/< 



{l + \x\^)^^+'^\D'u\'^dx 



1/2 



(1.11) 



where \D^u\ denotes the maximum modulus of any partial derivative of order i of the 
function u. If the quantity (1.11) is finite then u is said to belong to the weighted Sobolev 
space H^. Assume for asymptotic flatness that gab — ^ab G Hg and that kab £ 
for s sufficiently large and —3/2 < 5 < —1/2. If a spacetime is asymptotically flat then 
Theorem 1.1 can be sharpened to say that there exists a local solution corresponding to the 
given initial data and coordinates deflned for that solution such that the solution exists on 
X [0, T) for some T > and the data induced on the hypersurfaces t =const. satisfy the 
same type of asymptotic flatness conditions as the initial data on the hypersurface t — 0. 
More precisely, there exist coordinates so that the data induced by the solution on each 
hypersurface of constant time belongs to the same weighted Sobolev space as the initial 
data. Moreover the restrictions of the functions goo + 1 and goa of the metric components 
to any hypersurface of constant time belong to weighted Sobolev spaces. A proof of this 
is sketched in Section 4. 

The local existence theorem can be supplemented by a statement that the solution 
depends continuously on the data. This is not stated in full generality here; only some 
statements for asymptotically flat data are given. Given a family of initial data which is 
bounded in a weighted Sobolev space with s sufficiently large and is such that the metric 
is uniformly positive deflnite, the coordinates above can be chosen so that the solutions 
corresponding to all data in the family exist on the same time interval [0, T), the weighted 
Sobolev norm of the data induced by these solutions on a hypersurface of constant time 
is uniformly bounded and the induced metric on one of these hypersurfaces is uniformly 
positive definite. Moreover, the restrictions of ^roo + 1 and goa to each hypersurface of 
constant time arc uniformly bounded in a weighted Sobolev space for all data in the 
family and goo is uniformly bounded away from zero. See Section 4. 

Suppose now that an asymptotically flat initial data set admits a group G of symme- 
tries, i.e. that a Lie group G acts on the manifold S in such a way that gab, kab and fo 
are preserved. Then the spacetime in Theorem 1.1 can be chosen so that it admits G as 
a symmetry group. More precisely, there exists an action of G on M by isometrics which 
preserves / and which restricts to the original action on S. To prove this, consider the 
geodesic 'j{p) through a point p & S orthogonal to in a spacetime with the given initial 
data. Let t(j)) denote the largest number such that, when 7(p) is parametrized by proper 
time, with p corresponding to the parameter value zero, this geodesic is defined on the 
interval {—t{p),t{p)). By what has been said above, the spacetime can be chosen so that 
there exists a number T > which is smaller than t{p) for each p & S. Since the spacetime 
is globally hyperbolic, each point can be joined to the initial hypersurface by a timelike 
geodesic of maximal length. If T is chosen sufficiently small then the geodesic is unique. 
A new spacetime can be defined as the open subset of the original spacetime where this 
distance is less than T. It will now be shown that the action of G on 5 extends to an 
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action on this new spacetime, whose underlying manifold will be denoted by M. Given 
a point g G M, let p be the point where the unique geodesic 7 of maximal length from q 
to the initial hypersurface S meets 5". Let (p : G x S ^ S denote the action of G on 5". 
For g E G and q G M, let 4>g{q) be the point which lies the same distance to the future 
of (j)g{p) along the future-directed geodesic starting orthogonal to S at (j)g{p) as q lies to 
the future of p along 7(p). This defines a mapping : G x M — > M by ^{g.,q) = 4>g{Q)- 
This mapping cf) is an action of G on M which restricts to 0. By the uniqueness part of 
Theorem 1.1, it must preserve Qa/s and /. 

The global theorems to be proved later make use of the concept of maximal hypersur- 
faces. A spacelike hypersurface in a spacetime is called maximal if its mean curvature tik 
is zero. If an initial data set is given which is maximal in this sense and asymptotically flat 
it is of interest to know whether the corresponding local solution of the Einstein- Vlasov 
system, whose existence is guaranteed by the above theorem, can be foliated by maximal 
hypersurfaces in a neighbourhood of the initial hypersurface. Given what has already been 
said about the local existence of asymptotically flat spacetimes, this can be proved using 
the implicit function theorem (cf. [16]). The time coordinate t above can be chosen so 
that its level hypersurfaces are maximal hypersurfaces and it can also be arranged that as 
|a;| — > 00 the coordinate t agrees asymptotically with proper time along a geodesic which 
starts normal to the initial hypersurface. (This implies that qoq tends to unity as tends 
to inflnity.) When this has been imposed the foliation by maximal hypersurfaces is unique 
and so this construction gives a unique preferred time coordinate. The restriction of the 
solution to any of the maximal hypersurfaces is asymptotically flat. Once again there is 
a statement of uniformity. For a family of maximal initial data which is bounded in a 
suitable weighted Sobolev space the maximal foliation can be assumed to exist on a time 
interval which is uniform for all data in the family. These statements about the existence 
of maximal foliations are not too dependent on the fact that the matter is described by the 
Vlasov equation. The essential property which is needed for the existence and uniqueness 
theorems is the strong energy condition. 

The approach to studying the global structure of asymptotically flat solutions of the 
Einstein-Vlasov system presented in these lectures is closely related to work which has 
been done in the spatially compact case [25, 4]. The main difference is that in a non- 
flat spacetime satisfying the strong energy condition with a compact Cauchy hypersurface 
there exists at most one maximal hypersurface [16]. In this case the maximal foliation of 
the asymptotically flat case can be replaced by a constant mean curvature (CMC) foliation. 
This means that each leaf of the foliation has constant mean curvature, while the value of 
the mean curvature varies monotonically from one leaf to the next. 

2. Spherical symmetry 

Investigating the global properties of general solutions of the Einstein-Vlasov system 
is beyond the scope of existing mathematical techniques. For comparison, note that the 
same comment applies to the Vlasov-Maxwell system (cf. [19] for the most general known 
results) while general global existence results have been obtained for the simpler Vlasov- 
Poisson system ([17], [13], [21]). When a system of partial differential equations appears 
inaccessible to direct attack, a natural strategy is to study the simpler equations obtained 
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by imposing symmetry conditions on the solutions of the original equations. In the case 
of asymptotically flat solutions of the Einstein- Vlasov system, it seems that there are only 
two possible symmetry assumptions: spherical and rotational symmetry. In the latter 
case, where the symmetry group is one- dimensional and has fixed points (on the axis of 
rotation), the simplification obtained is not sufficient to bring the problem within range of 
present techniques. Thus in the following treatment of global questions we consider only 
the spherically symmetric case. Note for comparison that in the spatially compact case a 
wider variety of tractable symmetry types exists. 

A solution (M, gfQ,^, /) of the Einstein- Vlasov system is said to be spherically sym- 
metric if there exists an action of SO{3) on M by isometrics whose generic orbits are 
two-dimensional such that the natural lift of this action to the mass shell preserves /. 
There is of course an analogous definition of a spherically symmetric initial data set. Con- 
sider now a spherically symmetric asymptotically flat maximal initial data set. From the 
last section we know that there exists a corresponding local solution of the Einstein- Vlasov 
system. Moreover, it can be assumed without loss of generality that S0{3) acts on this 
local solution as a symmetry group so that it is spherically symmetric. Furthermore, there 
exists a neighbourhood U of the initial hypersurface which can be foliated by maximal 
hypersurfaces whose intrinsic geometry is asymptotically flat and this foliation is unique. 
It follows from the latter fact that each maximal hypersurface is invariant under the action 
of SO{3). In other words, it is a union of orbits of the action of SO{3) on spacetime. Let 
r, the area radius, be defined by the condition that on any orbit it takes the constant value 
A/4n, where A is the area of the given orbit. Consider now a fixed spacelike hypersur- 
face which is invariant under the action of SO{3). A geodesic of the induced metric on 
which starts orthogonal to the orbits remains orthogonal to them. These geodesies will 
be called radial geodesies. Let 9, be standard spherical coordinates on one of the orbits. 
Extend them to be constant along the radial geodesies. Since radial geodesies can never 
intersect except at the centre this prescription is globally well defined. In the coordinates 
(r, 9, 0) the metric components gi2 and gis vanish. The metric intrinsic to the orbits takes 
the standard form r'^(d9'^ + sin^ 9d4>^). Hence, provided the gradient of r does not vanish 
anywhere, the induced metric on Sq can be written in the form 



This way of defining spatial coordinates was used in [20] . It cannot be used if the restriction 
of r to the hypersurfaces of constant time has a vanishing gradient somewhere. A point 
where the gradient vanishes corresponds to a minimal surface. Here we use another type of 
coordinate system which does not suffer from this difficulty. First note that it is possible 
to write the metric in the form 



without restriction, for some function B. Since there is always a neighbourhood of the 
origin without minimal surfaces it can be assumed without loss of generality that r{x) = x 
near the origin. Since B is smooth when considered as a function on spacetime, it must 
be a smooth function of x^. Furthermore, in order that the spacetime be regular at a; = 



(2.1) 



B^x)dx^ + r^x)(d9^ + sin^ 9d(f)^) 



(2.2) 
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and not have a conical singularity, -B(O) must be equal to unity. This means in particular 
that it is possible to write B{x) = 1 + x'^D{x), where D{x) is a smooth fimction of x^. A 
new coordinate R will be sought which is a function of x and has the property that, when 
expressed in terms of the coordinate R, the metric takes the form: 

A^{R){dR^ + R^{de'^ + sin^ (2.3) 

Coordinates of this sort, or more precisely the Cartesian coordinates corresponding to these 
polar coordinates, are often known as isotropic coordinates. Writing down the coordinate 
transformation shows that R{x) is a solution of the equation dR/dx = RB{x)/r{x). This 
equation has a solution which is unique up to a constant scaling. To see this, note first that 
this is a first order homogeneous linear ordinary differential equation. Hence it has a one- 
parameter family of solutions which can all be got by multiplying one particular solution 
by an arbitrary constant. If we know the existence of the solution near the origin then 
the existence for all values of R follows by the standard existence and uniqueness theorem 
for ordinary differential equations. Near the origin the equation can be solved, giving 
R = Cxexp x' D{x')dx' . Near infinity, the asymptotic fiatness of the metric implies 
that A tends to a constant value. The scaling can be chosen so that lim^j^oo ^(-R) = 1 
and then the solution is determined uniquely. Now let t be a time coordinate which is 
constant on each leaf of the preferred foliation by maximal hypersurfaces and which agrees 
asymptotically with proper time. Introducing a coordinate R as above on each leaf puts 
the metric in the form: 

-a^ {t, R)dt^ + A^{t, R)[{dR + (3{t, R)dtf + R'^ide^ + siii^ ed(j)'^)] (2.4) 

As a consequence of the choice of coordinates the functions A and a tend to unity as 
R ^ oo for each fixed t while (3^0. The smoothness of the spacetime metric together 
with spherical symmetry implies that a, A, and R~^(3 are smooth functions of R^. 

For a metric of the form (2.4) which satisfies the extra condition that the hypersurfaces 
of constant time are maximal hypersurfaces the field equations and coordinate conditions 



take the following form: 

(i?2(Al/2)')/ = _1^5/2^2(|^2 ^ ^Q^p^ ^2.5) 

a" + 2a' /R + A-^A'a' = aA^^K^ + 47r(p + irS)] (2.6) 

K' + ?,{A-^A' + l/R)K = SttAJ (2.7) 

/?' - R-^f3 = laK (2.8) 

dtA = -aKA + (/3A)' (2.9) 
dtK = -A-^a" + A-^A'a' + a[-2A-^A" + 2A-'^A'^ - 2A-^A'/R - SttSr 

+ 47rtr5 - 47rp] + PK' {2 AO) 



The notation used in these equations will now be explained. A prime denotes a derivative 
with respect to R. The quantity K is that obtained by contracting the second fundamental 
form of the hypersurface t =const. twice with the unit vector A~^d/dR while Sr is 
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obtained in the corresponding way from the energy-momentum tensor. The quantity ItS 

is the trace of the spatial part of the energy-momentum tensor, i.e. if n'^ is the unit future- 
pointing normal vector to the hypersurfaces of constant time, then trS = T(^i3{g°'l^ + n"'n^). 
The quantity j is obtained by contracting Tap once with and once with the vector 
A~^d/dR and p is the energy density Tapn'^n^. In the standard terminology of the 3-1-1- 
decomposition of Einstein's equations, a is the lapse function and /3 is the one non-trivial 
component of the shift vector. Equation (2.5) is the explicit form of the Hamiltonian 
constraint (1.7) in this class of spacetimes with this kind of coordinate condition. The 
maximal slicing condition is expressed by the lapse equation (2.6). The one non-trivial 
component of the momentum constraint (1.8) in these spacetimes is (2.7). Equation (2.8) 
is a consequence of the coordinate condition chosen while (2.9) follows from the definition 
of the second fundamental form. Finally, (2.10) is the one non-trivial Einstein evolution 
equation in this class of spacetimes. This form of the field equations has been used by 
Shapiro and Teukolsky for numerical calculations (see [27]). 

What has been shown above implies that given asymptotically flat spherically symmet- 
ric maximal initial data for the Einstein- Vlasov system, there exists a corresponding local 
solution and some Ti > such that this spacetime can be covered by coordinates which 
cast it in the form (2.4) and for which the time coordinate ranges in the interval (— Ti, Ti) 
and the initial hypersurface is given by t = 0. In fact we are only interested in evolution to 
the future and hence only consider the part of spacetime on the interval [0, Ti). The quan- 
tities describing the metric and matter then satisfy the equations (2.5)-(2.10). The Vlasov 
equation is of course also satisfied. It turns out to be useful for some purposes to write 
it in Cartesian coordinates. Let be the coordinates (i? sin 6* cos 0, i? sin 6* sine/), i? cos ^). 
Define a related orthonormal frame by = A~^d / dx^ . Then if the mass shell is coordina- 
tized by (t, a;°',v*) , where denote the components of a vector in the given orthonormal 
frame the Vlasov equation takes the explicit form: 



dl 
dt 



+ aA~ 



X 



-(3 ^ 



dl 
dx 



+ 



1^— aK (v — 



R 



3Vr 



X ^ 

R^ 



-aA '^A' (vvr - Ivl'^T^') , ^ 

V ^ ^ R) yiTH 



dv 

(2.11) 

Here a dot denotes the usual inner product in R^, \v\ = ^Jv ■ v and = (v ■ x)/R. For 
given initial data there exists a solution on an interval [0, Ti) of the equations (2.5)-(2.11) 
supplemented by the definitions of the matter quantities. The whole system of equations 
will be referred to as the 'reduced Einstein- Vlasov system'. This solution of the reduced 
system is uniquely determined by its restriction to t = 0, as follows easily from the general 
uniqueness theorem for solutions of the Einstein- Vlasov system and the uniqueness of the 
coordinate system used. As a consequence there exists a greatest value of Ti (finite or 
infinite, call it T*) for which a solution of the reduced equations with the given initial data 
exists on the time interval [0,Ti). The interval [0,T^) is called the maximal interval of 
existence. The global existence question, which is the main theme of these lectures, is the 
question under what circumstances = oo. 

One possible strategy for proving global existence theorems will now be outlined. 
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Suppose that in some way it were possible to show for given initial data that for any- 
corresponding solution on a finite interval [0,Ti) the metric components, the distribution 
function and all their derivatives of all orders with respect to t and x were bounded. 
Then global existence for these initial data would follow. For the metric components, the 
distribution function and all their derivatives of all orders would be uniformly continuous. 
By a standard theorem on metric spaces they would all extend to continuous functions on 
the closed interval [0,Ti]. By another standard theorem, this time of real analysis, the 
extensions are and each derivative of each extension is equal to the extension of the 
corresponding derivative. In this way smooth initial data are defined on the hypersurface 
t = Ti. Provided these new initial data are asymptotically flat, the local existence theorem 
can be applied again to show that the original solution has an extension to an interval 
[0,T2) with T2 > Ti. Hence Ti 7^ T*. But since Ti was an arbitrary positive number this 
only leaves the possibility that = 00 and global existence is proved. In the following a 
situation is exhibited where bounds similar to those which are assumed in this argument 
can actually be obtained. 

To get closer to the situation which has just been described, consider a solution of 
the reduced equations defined on some interval [0,Ti). It will now be shown that many 
quantities can be bounded by using the Einstein- Vlasov system. By the dominant energy 
condition p > 0. Hence equation (2.5) shows that is a non-increasing function of 

R for each fixed t. When = it is zero and hence i?2(^V2y < q. It follows that A' < 0. 
The boundary condition that ^ — > 1 as — > 00 then gives ^ > 1. Next, from (2.6), 

(R^AaJ = aA^R'ilK^ + 4n{p + trS)] > (2.12) 

The expression on the right hand side of (2.12) is non-negative, as follows from the domi- 
nant energy and non-negative pressures conditions. Since R^Aa' vanishes for i? = it can 
be seen that R'^Aa' > and a' > 0. Using the boundary condition that o; — > 1 as — > 00 
gives a < 1. 

Next an estimate of Malec and O Murchadha [15] will be used. The expansions of the 
null geodesies which start normal to the orbits are given by 

d = 2{A-'A' + {AR)-') + K 

e' = 2(A-'^A' + {AR)-^)-K ^' ^ 

The area radius is given by r = AR. The theorem of [15] states that, if the dominant energy 
condition holds, the quantities r$ and r$' are bounded in modulus by two. Adding and 
subtracting these estimates and using the explicit expressions in (2.13) gives the inequalities 
< 2 and \RA~^A' + 1| < 1. Since ^ > 1 it can be concluded from the first of these 
inequalities that \K\ < 2R~^. The second inequality gives < 2R~^. In particular 

this gives pointwise bounds for K and A~^A' away from the centre. Integrating (2.8) gives 

R^^/3{R2) - Ri^/3{Ri) = ^ {aK/s)ds (2.14) 

Asymptotic flatness implies that K = 0{R~^) as i? — > 00. Hence it is possible to let R2 
tend to infinity in this equation to get an expression for R~^(3{R) as an integral from R 
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to OO. Using the bounds already obtained for K and a shows that 

/■OO 

\(3\ <lR {2/s^)ds < 3 (2.15) 
Jr 

and so is bounded. A bound for a' can be obtained by analysing the inequality 

{R'^Aa'y > 0, which was already used to show that a < 1. Integrating this between 
the radii Ri and i?2 with Ri < i?2 gives: 

a'{R2) > {Ri/R2)\A{R^)/A{R2))a\Ri) (2.16) 

Integrating again then gives: 

a(i?2) > a'(i?i)(i?f^(i?i)) f \-^{A{R))-^dR (2.17) 

JRi 

Now use the facts that a{R2) < 1 and that A{R) < ^(.Ri) for R> Ri to see that: 

rR2 

1 > a'{Ri)Rl I R-^dR (2.18) 

JRi 

This holds for all R2 > Ri and so it is permissible to replace the upper limit in the integral 
by infinity. Evaluating the integral gives a' (R) < R-^ for all > 0. Note that in deriving 

all these estimates, the only properties of the matter fields used were the dominant energy 
condition and the inequality p + trS > 0. The latter follows from the strong energy 
condition. Thus all these estimates hold not only for the Einstein- Vlasov system, but also 
for the Einstein equations coupled to any matter model which satisfies the dominant and 
strong energy conditions. This includes perfect fiuids with reasonable equations of state, 
matter described by the Boltzmann equation, the massless scalar field (or more generally 
wave maps), any of these matter models combined with an electromagnetic field in such a 
way that the total energy-momentum tensor is the sum of the individual energy-momentum 
tensors, and the Yang-Mills equations for any semi-simple gauge group. 

These estimates give good information on the solution away from the centre but, 
except in the case of the estimate for (3, give no control at the centre. Pointwise estimates 
for a'. A' and K which also give useful information at the centre can be obtained by 
optimization arguments. For any fixed radius Rq we have 

\K{R)\<A7rA\0)\\p\\^Ro (2.19) 

for any R < Rq. Here || ||oo denotes the L°° norm in space, i.e. the maximum value 
of a function on a hypersurface of constant time. On the other hand, if R > Rq then 
\K{R)\ < 2Rq^. Thus for any value of R it is true that 

\K{R)\ < 47rA\0)\\p\\^Ro + 2R-^ (2.20) 
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Thus 



This can be optimized by choosing Rq so that the function of Rq occurring on the right 
hand side of this last inequahty has a critical point. This occurs when Rq is equal to 

[27rA4(0)||p||oo]"^/^- It follows that ||i^||oo < CA2(0)||p||^^ for some constant C. A similar 
procedure can be used to estimate a'. 

a'{R) = A-^R-'^ + 47rp + 47rtr5)ds 

io (2.21) 

<CA^(0)||p||ooi?0 

for R < Rq. Combining this with the previous pointwise estimate for a' gives 

a' <C{A''{0)\\p\\ooRo + Ro^) (2.22) 

Doing an optimization as above then leads to an estimate of the form a' < C^^/^(0)||p||to^- 
Prom the equation for A: 

A'{R) = \R-^A^/^ / A^/^s^{lK^ + 167rp)ds (2.23) 
Jo 

\A'\ < C^(0)[i^6(0)||p||ooi?o + 2i?o '] (2.24) 

Optimizing gives ||^'||oo < C'74'^(0)||p(0)||to^- Starting from the same equations we can 
also bound R~^K, R~^a' and R~^A' pointwise in terms of ||p||oo and A{0). In particular 
R~^K can be bounded by a constant multiple of ^^(0)||p||. Equations (2.6) and (2.7) then 
allow a" and K' to be bounded. Solving (2.5) for A" shows that it too can be bounded. 
Equation (2.14) and the bounds for fi and R~^K imply that R~^P can be bounded by a 

1 /2 

constant multiple of A^(0)||p||oo • Using equation (2.8) then gives a similar bound for 

In the proof of variants of these estimates discussed in the next section, the con- 
servation of the total (ADM) mass plays a role and it is convenient to say something 
about this conservation law at this point. It is a quite general property of asymptoti- 
cally flat spacetimes. It is particularly easy to see in the situation considered here, where 
outside a compact set the spacetime is vacuum and spherically symmetric. Equation 
(2.7) can be rearranged to give {A'^R^K)' — SirR'^A'^j. In vacuum this integrates to give 
K = Ko{t)R~^A~^ for some function Ko{t). Putting this in (2.5) and using the vacuum 
condition again gives {R^{A^/'^yy = 0{R~^). This can be integrated to give: 

A{t, R) = {1 + Ao{t)R-^ f + 0{R-^) (2.25) 

It follows from (2.8) that (R-^/Sy = ^R'^aK = 0{R-^). Hence (3 = 0{R-^) and 
(3' = 0{R-^). It then follows from (2.9) that dfA = 0{R-^). Integrating this last 
equation in time from to t shows that A{t, R) — A(0, R) + 0{R~^), so that Ao{t) is in 
fact independent of t. The ADM mass is given by limji_,oo{—R^A') and so is also time 
independent. Prom (2.5) we can calculate that the ADM mass is equal to 

rriADM = I / A^/^R^i^K^ + 16TTp)dR (2.26) 
Jo 
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In all these estimates only the dominant and strong energy condtions have been used. 

To go beyond the results of the last paragraphs it is necessary to use the specific 
nature of the matter model. For the Einstein-Vlasov system a continuation criterion will 
be proved. It is formulated in terms of a quantity P{t), which is defined to be the largest 
momentum of any particle at time t. In other words 

P{t) = sup{\v\ : f{t, x,v)^0 for some x} (2.27) 

Before stating the continuation criterion it is necessary to take some time to discuss the 
relation between the differentiability of the functions of R describing the spacetime and the 
differentiability of the corresponding objects in spacetime. The simplest case is that of the 
scalar functions a and A. They are C°° in the spacetime sense if and only if they are C°° 
as functions of R and all the derivatives of odd order vanish at the origin. This follows from 
Lemma Al of the appendix with m — 0. There is also a quantitative version of this, which 
follows from Lemma A2. Consider next By definition (3 = (3°'Xa/R and = (3x°' /R, 
where (3°' is the shift vector. Because of spherical symmetry (3 must vanish at the origin. 
Hence we can apply Lemma A3 to Lemma A4 gives quantitative results for 

Consider next K. By definition A^K = kabx"'x^ /R"^ ■ The maximal hypersurface condition 
and spherical symmetry together imply that kab vanishes at the origin. Hence Lemma A2 
can be applied with m = 2. Similar considerations apply to the matter quantities j and 
Sr, whereby in the latter case it is necessary to write Sr = Sr + |tr5', with Sr being 
the contribution to Sr of the trace free part of Tab- Because of spherical symmetry Sr 
vanishes at R = 0. The following expressions for some of the quantities occurring in the 
Vlasov equation are also significant: 

pXa/R = Pa 
a'xa/R = VaO: 

kab = -^KA^iSab - 3XaXb/R^) 
A'xa/R = VaA 

Theorem 2.1 If a solution of the reduced Einstein-Vlasov system on the interval [0, T) for 
some positive real number T is such that P{t) and A{t, 0) are bounded then the solution 
extends to an interval [0, Ti) with Ti > T. In particular, if the maximal interval of existence 
[0, T*) is finite then either P{t) or A{t, 0) is unbounded there. 

Proof Suppose that P{t) is bounded on the interval [0,T). Then the matter quantities 
p, trS" and Sr are bounded there. It has been shown above that this, together with the 
boundedness of A{t^ 0), implies that the quantities A, A' , A" , R~^A', a, a' , a" , R~^a' , K, 
K' , R~^K. j3, P', and R''^(3 are bounded. It remains to show that all higher spacetime 
derivatives of all these quantities are bounded. We have a bound for a and A and a 

bound for K and /3 when these are considered as functions of R. These imply a 
bound for A, a and a bound for and kab in the three dimensional sense, using the 
results of the appendix. It follows that a bound for all the coefficients of the Vlasov 
equation on the support of / is obtained. The equations obtained by differentiating the 
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(2.28) 



Vlasov equation with respect to x and v then give the boundedness of the first derivatives 
of / with respect to x and v. Using the definition of the energy-momentum tensor gives 
a bound for its Cartesian components. The results of the appendix then imply a 
bound for p, j and trS". 

Higher derivatives can now be bounded inductively. Assume that a solution of the 
reduced equations on a given time interval is such that the C^"*"^ norms of a and A and the 
Qk j^oj-ms of 13, /, p, J, IyS and Sr are bounded and that is also bounded. Note 
that it has already been shown that under the hypotheses of the theorem this statement 
holds for /e = 1 and this suffices to start the induction. Now consider the case of general 
k. It is convenient to rewrite some of the reduced equations in the following form: 



Applying Lemma A5 to (2.29) and (2.30) gives C'^+i bounds for (A^S)' and Aa' . Com- 
bining this with the information already available gives (7^^+^ bounds for A and a. In a 
similar way, (2.31) and Lemma A5 give a C^^^ bound for K. The quantity /3'(0) is already 
known to be bounded. Hence (2.32) and Lemma A6 imply a C^'^^ bound for (3. (In fact it 
implies a (7^+^ bound, but that is not relevant here.) Using the relationships between dif- 
ferentiability of functions of three variables and functions of one variable discussed above 
and the results of the appendix, it follows that the coefficients of the Vlasov equation 
are C^'^^. Hence the solution of the Vlasov equation is bounded in the C^'^^ norm. An 
immediate consequence is that the Cartesian components of the energy-momentum tensor 
are bounded in the C^'^^ norm. Finally, applying the results of the appendix again shows 
that the (7^+^ norms of p, j and tvS are bounded and this completes the inductive step. 

It was mentioned earlier that in order to prove global existence it would suffice to bound 
the derivatives of all orders of all quantities of interest with respect to t and R. Here 
only the derivatives with respect to R have been bounded and it turns out to be difficult 
to boimd the time derivatives of the lapse function directly. Fortunately it is enough, in 
the present context, to bound the spatial derivatives, as will now be shown. Let t„ be a 
sequence of times with <T for each n and lim^^oo tn = T. The initial data induced by 
the given solution on the hypersurfaces t = tn define a sequence of initial data which are 
bounded in the C°° topology. In fact they are also bounded in the topology of a weighted 
Sobolev space. To prove this it suffices to obtain some estimates on an exterior region, 
say that defined by i? > 1. Equation (2.29) and the conservation of ADM mass shows 
that A' is 0{R~^), uniformly in t. Equation (2.31) and the fact that the support of the 
matter is contained in a region of the form R < Rq for all t in the interval [0, Ti) shows that 




(2.29) 



(2.30) 



(2.31) 



(2.32) 
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K = 0{R~^), uniformly in t. Using (2.25) and the fact that under the given circumstances 
the 0{R~^) error term there is uniform in t shows that A — 1 = 0{R~^), uniformly in 
t. Thus Qab — Sab IS bouudcd in Hg and kab is bounded in Hg^-^ for —3/2 < d < —1/2. 
To apply the more precise version of the local existence theorem for asymptotically flat 
spacetimes it is necessary to have a similar statement for weighted Sobolev spaces of higher 
order. This can be proved straightforwardly by induction using the equations (2.5) and 
(2.7). It follows that the solutions of the Einstein- Vlasov system corresponding to the data 
on each of the hypersurfaces of constant t exist on some time interval of length e about 
the initial time where data are given, with e independent of n. Hence the solution extends 
to the interval [0, T+e). 

With this result in hand, it is easy to show that the first singularity, if one exists, must 
occur in the centre. 

Theorem 2.2 If a solution of the reduced Einstein- Vlasov system on the interval [0, T) for 
some positive real number T is such that it has a smooth extension to an open neighbour- 
hood of the point with coordinates (T, 0) then the solution extends to an interval [0,Ti) 
with Ti > T. In particular, if the maximal interval of existence [0, T*) is finite then the 
solution has a singularity at the point (T*, 0) 

Proof The neighbourhood occurring in the hypotheses of the theorem contains all points 
with t > Ti and R < Ri for some T\ <T and some J?i > 0. Since the solution is smooth 
for t < Ti it follows that all unknowns in the reduced system are bounded on the region 
R < Ri, < t < T. The value of A at any point of a hypersurface of constant time can be 
bounded by its value at the centre at the given time and so A is bounded on the interval 
[0, T). It was shown earlier that on any region of the form R > Ri the quantities A~^A' 
K and a' are uniformly bounded. Hence under the present assumptions the quantities A' , 
K and a' are bounded on the interval [0, Ti). It was also shown that /3 is always bounded 
everywhere. It can be concluded that the coefficients in the Vlasov equation are bounded 
everywhere. Hence P{t) is bounded on the interval [0,Ti), and applying Theorem 2.1 
completes the proof. 

Remark It is clear from the proof that the existence of an extension could be replaced 
by the assumption that there exists some i?i > such that p and A are bounded on the 
region R < Ri 

Theorems 2.1 and 2.2 are analogues of Theorem 3.2 of [20] and Theorem 4.1 of [22] 
respectively. It is instructive to compare the theorems involving maximal-isotropic coordi- 
nates proved here with those involving Scliwarzschild coordinates proved in [20] and [22]. 
The continuation criterion of Theorem 2.1 appears at first sight weaker than that of [20] 
since it is assumed that not only P{t) but also A{t, 0) is bounded. However, it is much 
easier to pass from Theorem 2.1 to Theorem 2.2 than it is to pass from the continuation 
criterion of [20] to the regularity theorem of [22]. Moreover, the passage from Theorem 2.1 
to Theorem 2.2 does not involve any deep analysis of the Vlasov equation, which the proof 
of Theorem 4.1 of [22] does. Thus it is reasonable to hope that the method of proof used 
here can more easily be adapted to matter models other than the Vlasov equation than 
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the approach using Schwarzschild coordinates. In the next section it will be seen that the 
apparently weaker continuation criterion given by Theorem 2.1 is also good enough to be 
applied in the proof of a global existence theorem for small initial data. 

At this point some further remarks on the notion of 'well-behaved' matter models are 
in order. Consider the case of dust, i.e. a perfect fluid without pressure. In fact (see 
[24]), smooth solutions of the Einstein-dust equations can be considered as distributional 
solutions of the Einstein- Vlasov system. Dust satisfies the dominant energy and non- 
negative pressures conditions. However it cannot be expected that an analogue of Theorem 
2.2 holds for dust. The reason is the occurrence of so-called shell-crossing singularities, 
which do not occur at the centre. As has been discussed in [24] and [23] this is of significance 
for the formulation of the cosmic censorship hypothesis and Theorem 2 can be taken as 
an indication that the Einstein-Vlasov system is a good starting point for studying the 
cosmic censorship hypothesis and has advantages over other, superficially simpler matter 
models, such as a perfect fiuid. This is one of the main motivations for investigating the 
global properties of solutions of these differential equations. 

3. Global existence for small data 

In the last section a continuation criterion was given for solutions of the reduced 
equations. Now it will be applied to obtain a global existence theorem in a particular 
situation, namely that of small data. The notion of smallness of initial data is defined in 

the present context in terms of three quantities which characterize the size of the data. 
Let Fo = 11/(0) II oo, Po = -P(O) and let Rq be the smallest value of R such that /(O, i?) 
vanishes for R> Rq. 

Theorem 3.1 Let K be a fixed positive constant and consider initial data for the reduced 
equations with Rq < K and Pq < K. Then there exists an e > such that for all data of 
this type which, in addition, satisfy Fq < e the corresponding solution exists globally in 
time and the spacetime which it defines is timelike and null geodesically complete. 
Remarks 1. The spacetimes of the theorem are also spacelike geodesically complete but 
this will not be proved here. It is the completeness of timelike and null geodesies which is 
most interesting physically, since these represent the wordlines of particles. 
2. The statement on geodesic completeness is an important part of the theorem since 
a theorem on global existence in some coordinate time does not necessarily imply any 
interesting invariant information. 

In fact more detailed information concerning the asymptotic behaviour of the spacetimes 
covered by the theorem will be obtained. In particular, information will be obtained on 
the decay of the curvature as t — > oo. A good understanding of the decay properties 
of the curvature is also important for the proof of the theorem and for this reason the 
curvature components in a Cartesian frame will now be examined in some detail. The cur- 
vature can be computed using the following relations (which are independent of symmetry 
assumptions): 

Rabcd — Pabcd + ^ackhd ~ kad^bc 
^^^Rabcdn'' = -Vckab + Vbkac (3-1) 

Raarbn" = -87r[S'a6 + |(p - irS)gab] + Rab + trA;A;a6 - KcK 
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The first of these equations is the Gauss equation, the second the Codazzi equation and 
the third the Einstein evolution equation. It turns out that the frame components of the 
curvature tensor which are purely spatial are linear combinations of the quantities A~^A", 
^-4^/2 j^2^ with coefficients which are homogeneous functions of the Cartesian coor- 
dinates. The frame components with two indices corresponding to the time direction are 
linear combinations of these quantities and frame components of the energy-momentum 
tensor. The frame components with precisely one index corresponding to the time direction 
need to be calculated explicitly. They are given by: 

Roijk = lR-\xj5ik - Xk5ij){K' + ZKR-^ + ZA'^A'K) (3.2) 

The combination of metric coefficients which occurs is precisely that which is familiar from 
the momentum constraint. 

The idea behind the proof is as follows. In flat space free particles which start in a 
compact set spread out linearly with time. This causes the associated density to decay. 
In fact it decays uniformly in space like as t — oo. It is reasonable to suppose that 
a spacetime which evolves from 'small data' (i.e. data which are 'close' to data for flat 
space) has small curvature, so that the behaviour of solutions of the Vlasov equation is 
similar to that in flat space. Thus, with luck, the density will fall off at the same rate as in 
flat space. Conversely, it is this fall-off of the density which ensures decay of the curvature. 
The proof which follows makes these intuitive considerations precise and quantitative. It is 
broken up into a number of lemmas which are arranged in a way which parallels as closely 
as possible the proof of the analogous theorem in [20] . 

Lemma 3.1 Consider a spherically symmetric solution of the reduced equations on a 
time interval [0, T) with the property that the support of the restriction of T'*^ to each 
hyper surface t =const. is contained in the ball of radius Rq -\- 1 about the centre and 
^(^5 0) < 3 and suppose that: 

||p|U(l+t)'+'<i^i (3.3) 

for some constants Ki > and 5 G (0, 1]. Then there exists a constant C, only depending 
on Ki, S and the restriction of the solution to the initial hypersurface, such that: 

\\A"\\^{1 + tf+' + P'|U(1 + t)^+' + ||a'||oo(l + t)^+' + \\K\\oo{l + ty+^ < C (3.4) 

Moreover the L°° norm of the frame components of the curvature tensor can be bounded 

by a constant times (1 -|- t)~'^~^ . 

Proof Equation (2.29) and the conservation of the ADM mass (2.20) imply that A'{R) 
can be bounded by an expression of the form CRq^ for some constant C for all R > Rq. 
The constant depends a priori on a pointwise bound for A{t, 0) but this is taken care of 
by the assumption on A(t, 0) occurring in the hypotheses of the lemma. An analogous 
estimate holds for a' , as a consequence of (2.30). In the case of K it is necessary, in order 
to get an estimate which holds inside the matter, to cancel one power of R with a power 
of s in (2.31). Thus once again the estimate obtained involves Rq^- In the last section 
estimates in terms of Rq^ were used in an optimization argument to bound A', a' and K 

by IIpIIoo • If they are replaced by the estimates in terms of Rq just discussed, then the 
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optimization argument allows the L°° norms of A' , a' and K to be bounded by a constant 

2 /3 

multiple of ||p||oo • This is enough to take care of the last three terms in (3.4). It was 
remarked in the last section that R~^A' can be bounded pointwise by a constant times 
IIpIIoo if ^(^, 0) is known to be bounded. Putting this into equation (2.5) together with 
the information just obtained gives the desired estimate for the first term in (3.4). The 
estimate for the components of the curvature tensor follows from the expressions for these 
components given above. 

Remark The estimates obtained by optimization in the proof of this lemma could have 
been used in Section 2 instead of the other estimates obtained by the same method which 
were actually used there. The reason for presenting both types of estimates is that while 
the estimates in the proof of the lemma are stronger in situations where the solutions are 
small, the estimates of Section 2 are stronger where the solutions are large. Thus they 
might be important in other contexts. 

The following simple lemma is taken directly from [20] . 
Lemma 3.2 Consider the ordinary differential equation du/dt = F{t, u) for a function 
F satisfying the inequality \F{t, u)\ < r]{l + t)~^~^{l + \u\) for some constants 5 > and 
77 > 0. Then given any initial datum at t = the corresponding solution exists on the 
whole of [0, oo) and satisfies the inequality \u{t) — tt(0)| < | exp ^(1 + |w(0)|). 

This lemma is used to control the behaviour of timelike geodesies in a spacetime satis- 
fying certain inequalities. Consider a timelike geodesic 7 passing through a point with 
coordinates (t, x") and suppose that it intersects the initial hypersurface t = 0. Let r be 
proper time measured along 7, starting at t = 0. Let {e^} be an orthonormal frame which 
is parallelly transported and is such that Cg is the tangent vector to 7. Let 9''^ be the 
dual coframe. A Jacobi field along 7 is the derivative with respect to the parameter of a 
one-parameter family of geodesies in which 7 is embedded. If it is expressed as a linear 
combination Z^e'g of e'^, e'2 and 63 then satisfies the equation (see [12], p. 96): 

d'zydT^ = (i^^^^Ce^e^e^f )Z* (3.5) 

It turns out to be crucial for the global existence theorem to estimate the Jacobian deter- 
minant of the mapping 

^X"(0,t,a;«,^;^) (3.6) 

for fixed values of t and x", where X"'{s, t, x"', f *) is part of the solution of the characteristic 
system, as discussed in Section 1. Note that, as indicated by the use of the notation v^, 
the characteristic system which is of interest here is that of the Vlasov equation written in 
terms of frame components. This mapping can be described in words as follows. Follow 
the geodesic 7 through (t, x"') with initial tangent vector v^e^ backwards until it meets 
the hypersurface t = 0. The derivative of this mapping takes vectors tangent to the mass 
shell at the point with coordinates to vectors tangent to the initial hypersurface. A 
vector tangent to P at c'q can be identified with a vector in the tangent space to M at 
(t, x"') which is orthogonal to Cg. This can then be expressed as a linear combination Y^e'g. 
The frame vectors can be fixed uniquely by requiring that they be obtained from the 
coordinate vectors d/dv'^ by means of the Gram-Schmidt process. The Jacobian of the 
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mapping which transforms from the basis djdv^ to the basis and that of its inverse 
can be bounded in terms of \v\. The derivative of the mapping (3.6) can conveniently be 
written as a composition of three Hnear mappings from to itself. The first mapping 
Li is the one sending the components of a vector tangent to the mass shell in the basis 
{djdv"-^ to the components . Let L2(F^) be the value at t = of the solution of 
(3.5) with the initial data = and dZ^jdr = at the point {t,x''). Let LsiZ^) be 
the Cartesian components of the vector obtained by projecting the vector Z^e'g onto the 
hypersurface t — along the vector Cq. The components of this vector are given explicitly 
by: 

W^« = Z^[e--(ef/e[,")e-] (3.7) 
the derivative of (3.6) at the point of interest is L3L2L1. 

Lemma 3.3 Consider spherically symmetric solutions of the Einstein- Vlasov system on 
intervals [0, T) satisfying the following conditions: 

Il/olloo < e, 

Po < i^2, } (3.8) 



Ro < Ks 



the inequality (3.3) and the condition that ^ < 3. Then if e and Ki are sufficiently small 
there exists a constant K4, depending only on e, Ki, K2 and ^3, such that a > 1/2 and: 

\\Ta(3\\oo<K4{l+t)-^ (3.9) 

Proof Lemma 3.1 provides a lot of information on the decay of geometric quantities as 
t 00. In particular the estimate (3.4) can be combined with Lemma 3.2 to show that if 
Ki is small enough 

Pit) < Po + 1 (3.10) 

It follows that a'{R) < Ce{R + R'^). If e is small enough then this can be integrated to 
show that a > 1 — 3Ce/2. This implies the first statement of the theorem. The inequality 
(3.10) implies a uniform bound for the Cartesian components of Cq if the tangent vector 
to the geodesic 7 is contained in the support of the distribution function. (Only geodesies 
of this kind are of interest here.) The remaining frame vectors satisfy: 

+ r^^ef e;7 = (3.11) 

The spatial components of (3.11) can be written as 

d 



dr 



{e'n = -r^-yef 6(7 (3.12) 



To profit from this it is necessary to have some information about the relation between 
proper time and coordinate time along 7. In fact 

dT/dt = a{l + \v\^)-^/^ (3.13) 
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Since a < 1 the relations (3.12) and (3.13) imply that 



<C(|e'/| + K°|)(l + t) 



-1-5 



(3.14) 



On the other hand, since is a unit vector, le^'^l < C|eg"|. Consider now the vectors 
along all geodesies contained in the support of /. We have a uniform bound for the initial 
data for equation (3.12) and so (3.14) and Lemma 3.2 imply the global boundedness of e'^ 
and hence of e'g. (Global boundedness means by definition that they can be bounded by 
a constant which is independent of T.) 

Consider now once again a timelike geodesic 7 whose tangent vector is contained in 
the support of the distribution function. As a consequence of (3.13) there are positive 
constants Ci, C2 such that 

Cit<T< C2t 



along 7. We can assume without loss of generality that C*2 > 1 and then 
Thus if — R^^^O'^Cq Cq ^ an estimate of the form 



|K/(r)|<C(l + r) 
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(3.15) 



(3.16) 



(3.17) 



holds along 7. Let tq be the value of r at the point (t,a;") and let Z^{t) be the solution 
of (3.5) with Z'{to) = and rfZVdr(ro) = Y\ Let 



E\t) = Z%t)-{t-to)Y^ 
Then by Taylor's theorem with integral remainder 

E\t)= I {To-a){a-T)K^{a)Y*da+ f {t - a)K^{a)E\a)da 

J To J To 

The first integral can be estimated as follows: 



(3.18) 



(3.19) 



^ \T-a){a-To)Kt{a)Y'da < {tq - t)\Y'\ jj Ca{l + a 



-2-5 



da 



(3.20) 



<C(To-T)|r* 



since a{l + cr) ^ ^da < 00. Combining (3.19) and (3.20) and applying Gronwall's 
inequality gives 



To 



Hence 



|£^'(t)| <C(ro-r)|F*|exp j C(a-r)(l + a) 

\E^{t)\<C{to-t)\Y'\ 
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da 



(3.21) 
(3.22) 



If Ki is chosen small enough then C < ^ and using (3.22) and the definition of shows 
that I dct £2! > Ctq. As a consequence of (3.15) Tq can be replaced by t in this inequality. 
The determinant of Li can be bounded from below uniformly for all geodesies whose 
tangent vectors are contained in the support of the distribution function. The same is true 
of the determinant of L3, as can be seen from (3.7). Thus the statement is obtained that 

> Ct^ (3.23) 

This has been shown for all v"" in the support of the distribution function. If we knew that 
the mapping (3.6) was injective on the support of / then we could change variables from 
v'^ to X"- in the definition of the energy-momentum tensor. As a consequence of (3.23) and 
the boundedness of P this would give an estimate of the form (3.9), completing the proof 
of the lemma. 

It will now be shown that for Ki sufficiently small this mapping is indeed injective. 
Suppose that, on the contrary, there are two distinct geodesies 70 and 71 starting at the 
point with coordinates (t, x") which meet the initial hypersurface at the same point. Let Vq 
and vl be the components of their initial tangent vectors on the mass shell. For / e [0, 1] let 
vl = {1 — 1)vq + lv\. Let 7/ be the geodesic with initial tangent vector corresponding to v]. 
Denote the spatial coordinates of the point of intersection of 7; with the initial hypersurface 
t = by ^"{l). Then ^«(0) = ^(1) and ^"(0 is a closed curve. Let = v\ - v'q. The 
tangent vector to the curve ^""{l) is the image under the mapping L3L2-L1 of w\ It can 
be concluded from (3.23) that the tangent vector to the curve ^"(/) can never vanish. It 
is convenient at this point to make a different choice of the frame e^. The change is that 
instead of the basis {d/dv'^} of the tangent space to the mass shell as above, a basis is 
chosen which consists of vectors which are linear combinations of the vectors djdv^ with 
constant coefficients with the first vector of the basis being w'^djdv^. The vector L\{vS) 
is proportional to (1,0,0). Using the smallness assumption on the data shows that the 
components e'f (0) are close to the components e'j^^(ro). More precisely, for K\ small enough 

|e;'^(0)-e;'^(ro)|<ry|e;«(To)| (3.24) 

for any given ry > 0. Similarly 

|e[r(0)-e[r(To)|<r7|e[r(To)| (3.25) 

In the same sense L2Li{w) is close to a vector of the form (^^,0, 0). For is close to 
and the frame vectors have also only changed a little between t — tq and r = 0. The 
tangent vector to the curve C"(0 given by the projection of L2Li{w) on the hypersurface 
t = along the vector eQ(0). The explicit form of this projection is displayed in equation 
(3.7). It will now be shown that, for Ki sufficiently small the second term in the square 
brackets in (3.7), with s = 1, can be bounded in modulus by a constant k < 1 times the 
modulus of the first term. Before proving this statement it will be shown that it implies the 
desired result. is as close as desired to a vector of the form [Z^ , 0, 0) and so W"- is as close 
as desired to the projection of a vector proportional to e'^ (0). The estimate for the terms in 



det 



(0,t,a;«,-) 
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(3.7) with s = 1 implies that the projection of a vector proportional to is contained in a 

convex cone about ei . Hence VF" is contained in a convex cone about a vector proportional 
to w°'. This is inconsistent with the fact that ^"(/) returns to its starting point. Thus it only 
remains to prove the above statement about the relative sizes of the terms in (3.7). Now 
e'l is close to A~^Bw°' for a certain quantity B, e'^ is close to Ba~^{w ■ vi){l + If/p)"-*^/^, 
s'q is close to A~^vf and is close to a~^(l + |fzp)^/^. Thus the two terms are close to 
Bw"" and B{w ■ vi){l + \vi\'^)~^v1 . Using the Cauchy-Schwarz inequality the modulus of 
the second expression can be bounded by |w| |t';p(l + I'WzP)"^ while the modulus of the 
first is \B\\w\. Using the fact that \vi\ < P, so that \vi\'^ / {1 + \vi\)^ < P^/{l + P^) <k <1 
completes the argument. 

Lemma 3.4 Suppose that initial data for the reduced system satisfy the inequalities (3.8) 
and let Tq be a fixed positive number. Then for fixed K2, there exists e > such that 
the solution corresponding to this data exists on the interval [0,To], and P{t) < 2P(0), 
^(^5 0) < 2 and ||p||oo is as small as desired there. 

Proof First a system of integral inequalities for P{t) and Q{t) = A{t, 0) will be derived. It 
is elementary that \\p\\oo < (^11/1100(1 + -P(^))^ with a constant C independent of the initial 
data. It was shown in the last section that the coefficients in the part of the characteristic 

1 /2 

system controlling the velocities can be bounded by C||p||oo (1 + Q(^))'^, for some positive 
integer q. Hence 

P{t) < P(0) + [' P{s){l + P{s)f{l + Q{s)yds (3.26) 

Equation (2.9) and the estimates obtained above imply that 

Q{t)<Q{0) + Ce^/^ [ P{s)il + Pis)Yil + Qis)yds (3.27) 

for some positive integer p. Let zi and Z2 be the unique solutions of the system of integral 
equations obtained by replacing the inequalities in (3.26) and (3.27) by equalities and P 
and Q by Zi and Z2 respectively, with initial data Zi{0) = -P(O) and -^2(0) = Q(0). When 
e is zero the solution of the integral equations is constant. In particular, it is global in 
time. Moreover the initial data zi{0) and 2:2(0) can be bounded in terms of i^2, K3 and 
e. It follows that for e sufficiently small the solution of these integral equations exists on a 
time interval [0, Ti) with Ti > Tq and, by making e smaller if necessary, it can be assumed 
that zi{t) < 2P(0) and Z2{t) < 2 for < t < Ti. Comparing the integral inequalities 
with the integral equations shows that P{t) < Zi{t) on any interval where both are defined 
and similarly Q{t) < Z2{t). By the above estimates ||p||oo will be as small as desired if e 
is chosen sufficiently small. The continuation criterion given by Theorem 2.1 implies that 
the solution can be extended to any time interval where the solution (zi, Z2) of the integral 
equations exists and this completes the proof of the lemma. 

Theorem 3.2 Let a non-negative C°° compactly supported spherically symmetric maxi- 
mal initial datum /o for the Einstein- Vlasov system be given which satisfies the inequalities 
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(3.8) for some positive constants e, K2 and K^. Then if e is small enough the corresponding 
solution of the reduced system exists globally in time. Moreover, for this solution 

\\To.(i\\oo<C{l+t)-^ (3.28) 

the metric coefficients and A are bounded and the estimates (3.4) hold with 5 = 1. 
Proof Let Ki and e be positive constants which are small enough so that the conclusions of 
Lemma 3.3 hold for some 6 < 1. Let Ti be a positive number satisfying K4^{l+TiY~^ < Ki. 
By Lemma 3.4 the constant e can be chosen so small that the solution exists on the time 
interval [0,Ti). Moreover, it can be arranged that on this interval the solution is as small 
as desired. In particular, e can be chosen so that (3.3) is satisfied on the interval [0, Ti) 
and A < 2 there. Consider a fixed initial datum satisfying (3.8). Define T^, to be the 
supremum of those positive numbers T such that the solution of the reduced equations 
corresponding to the given initial data exists on [0,T), and A < 2 for the solution on this 
interval and it satifies (3.3) there. Here = 00 is possible. In fact we will show that it is 
the only possibility. For suppose that < 00. The assumptions already made ensure that 
> Ti. The definitions of Ti and and Lemma 3.3 then show that the continuation 
criterion is satisfied on [0, T*) and that 

||T«^(T.)||oo < Ki(l + T,)-2-'5 (3.29) 

This means that the solution can be extended to an interval [0,T2) with T2 > T*. Also 
T2 can be chosen so that (3.3) is satisfied there. This contradicts the definition of T* and 
so in fact it must be the case that = 00. The inequality (3.3) holds on [0, 00) for some 
5 < 1. Applying Lemma 3.1 shows that an inequality of this form also holds for 5 = 1. 
The remaining conclusions of the theorem then follow from Lemma 3.1. 

The statement of Theorem 3.2 includes all the conclusions of Theorem 3.1 except that 
concerning geodesic completeness. Equation (3.15) shows that along a timelike geodesic 
proper time and coordinate time are equivalent. This, together with global existence in 
coordinate time shows that timelike geodesies are complete. Similarly, the information 
which we have on the geometry is more than enough to show that an affine parameter 
along a null geodesic is equivalent to coordinate time. 

At the end of the last section some remarks were made about the relation of the 
Vlasov equation with dust. For the Einstein-dust system there is no reasonable smallness 
assumption on initial data which ensures geodesic completeness of the corresponding so- 
lutions, as follows from the work of Christodoulou [5]. The only other matter model for 
which a global existence theorem for spherically symmetric asymptotically fiat solutions 
of the Einstein-matter equations with small initial data has been proved is the massless 
scalar field [6] . The method of proof used here in the case of the Vlasov equation can not 
be applied directly to the case of matter models with radiation such as the massless scalar 
field, since there the fiat space fall-off rates for the matter only make the energy-momentum 
tensor fall-off like and do not furnish the faster decay rates used in the above. 
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4. Local existence in general 

The purpose of this section is to present some aspects of the local existence theorem 

for the Einstein- Vlasov system, without any symmetry assumptions. The standard method 
for proving such theorems proceeds in several steps. First, some coordinate conditions are 
imposed, leading to a system of 'reduced equations'. This is similar to what was done 
in Section 2 in the spherically symmetric case. In the general case different coordinate 
conditions are used and hence the reduced system is also different. The second step is to 
prove a local existence theorem for the reduced equations. The third step is to establish the 
connection between the reduced equations and the full equations. This means showing that 
if the data for the reduced system satisfies certain gauge conditions, if the constraints are 
satisfied on the initial hypersurface and if the reduced equations are satisfied everywhere, 
then the coordinate condition and the constraints are satisfied everywhere. This then 
implies that the solution of the reduced equations is actually a solution of the Einstein 
equations. The first and third steps are not discussed further here; the reader is referred 
to [8] for details. (The treatment which follows uses the harmonic coordinate condition, 
which is that discussed in [8].) 

The reduced equations in harmonic coordinates constitute a system of nonlinear wave 
equations for the metric coupled to the Vlasov equation for the distribution function /. 
Solutions of this system have a domain of dependence determined by the light cone and so 
when proving a local existence theorem spatial boundary conditions play no role. For this 
reason it is sufficient to consider data of compact support on R"^. The data for the reduced 
system consists of all components of the metric and their first time derivatives together 
with the distribution function. It is assumed that gai3 — rjai3, dtg^p and / have compact 
support on the initial hypersurface. Such data can never satisfy the constraints globally 
on R"^ except in the case of data for flat space but this is not a problem, because of the 
possibility of using the domain of dependence to localize. To start with only the case of 
C°° data is treated but it will be shown later that it is easy to extend the argument to 
data of finite differentiability. 

Before going further some remarks will be made on the various notations used for 
derivatives in this section. This is intended to avoid confusion which might arise from 
the mixture of notation from differential geometry and from analysis which occurs. As in 
previous sections Greek or Roman indices attached to geometric objects are usually to be 
thought of as abstract indices although they may also occasionally denote components in 
a coordinate system. No confusion should arise from this dual role. Indices of this kind 
are also used to label the coordinate functions themselves and the partial derivatives with 
respect to the coordinates, denoted by d. The ranges of the indices are as before. On the 
other hand the indices on the operator D denoting differentiation are of a different kind. 
Here use is made of the multi-index notation which is very effective in handling expressions 
containing high order derivatives of functions of several variables. These multi-indices 
always refer to derivatives with respect to the spatial coordinates x". The expression D"u 
denotes a particular higher order partial derivative of the function u with repect to the 
spatial variables. The order of this derivative is denoted by \a\. Roman indices on D are 
used in connection with norms as in equation (1.11). We have 11-0*^11 = max|Q,|=j HD^mH. 
This notation may be applied to any norm. 
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In order to have a suitable framework for proving a local existence theorem it is 
necessary to work with function spaces which are well adapted to the equation being 
studied. The natural spaces for hyperbolic equations are the Sobolev spaces. In the present 
context, where the asymptotic behaviour plays no role, it is possible (and convenient) to 
use Sobolev spaces without weights. If tt is a smooth function on define: 

1/2 

(4.1) 

Compare this with the weighted Sobolev norms defined in (1.11). 
The reduced Einstein equations take the form: 

= F^p{g^s, d^g^s) + STrfT^,^ - ^^^T^sgap] (4.2) 

It is convenient (although not essential) to write this in first order form by introducing the 
first derivatives of the metric as additional variables. Let hafB-y = d^gap. Then in terms of 
the unknowns gap and hap-y the equations can be written as follows: 

—g^^dahapQ — 2g'^°'daha/3o = g°'^dahocpb + • • • 

^"'^O/Ia^a = g'^'dahapo (4-3) 

dogap — hapo 

The terms which are not written out explicitly are those coming from the right hand side 
of (4.2). The first of these can be written as a function of g^p and h^p-y without using 
derivatives. The reason why certain linear combinations of the original equations have been 
taken when writing (4.3) is to make contact with the notion of a symmetric hyperbolic 
system. Let tt be a function which takes values in an open set U of R'^ for some k. A 
differential equation for u of the form: 

A^{x", u)dou + A"(x", u)dau + B{x'^, u) = (4.4) 

is called symmetric hyperbolic if the matrices A^{x°',u) and A'^{x°',u) are symmetric for 
all (x'^ju) and if ^°(a;",tt) is positive definite. Here A^, A"' are functions on the product 
of U with some open subset of R"^ with values in the k x k matrices and S is a function 
on the same domain with values in R'^. li u — {gap — TjaPi hap-y), with rjap denoting the 
components of the Minkowski metric in standard coordinates then the equations (4.3) are 
of this form, for each fixed Tap. Moreover the reduction process, which is not described 
explicitly here, can be (and usually is) done in such a way that goo = —1 and goa = on 
the initial hypersurface. The open set U is defined by the condition that Qap have Lorentz 
signature. The Vlasov equation will be written in the form (1.1). Note that the Christoffel 
symbols are rational functions of the unknowns u. The full system of reduced equations 
consists of (4.3), (1.1) (with the Christoffel symbols expressed algebraically in terms of 
u) and the definition (1.3) of Tap- The Vlasov equation itself is a symmetric hyperbolic 
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equation of a particularly simple type, with only one unknown. This symmetric hyperbolic 
equation is defined on rather than R"^. The main theorem can now be stated: 
Theorem 4.1 Let {uq, /q) be C°° compactly supported initial data for the reduced 
Einstein- Vlasov system written in first order form. Then there exists a T > and a 
C°° solution (w, /) of the reduced system on the interval [0, T) which induces the given 
initial data. Moreover f{t) has compact support for each fixed t. 

Remarks 1. The domains of definition of the functions tto, /o, u and / are R^ , R^ , 

R3 X [0, T) and R^ x [0, T) respectively. 

2. The solution {u, f) is the unique solution with the properties stated in the theorem. 

In proving this result an existence theorem for solutions of linear symmetric hyper- 
bolic systems with C°° initial data is assumed (see e.g. [9], p. 668). In order to prove 
Theorem 4.1 a certain iteration is defined and then it is proved that this iteration converges 
in an appropriate sense. The convergence follows from certain inequalities satisfied by the 
solutions of linear symmetric hyperbolic equations. The proof makes use of the following 
Moser-type inequalities [14, 1] . Here these inequalities are only needed in the case of func- 
tions which are smooth and compactly supported. The first inequality concerns products 
of functions and says that if / amd g are smooth functions of compact support on R"^ then 

\\D'{f9)h < CdpVlblblloo + ll/lloolP^^lb) (4.5) 

Here || ||2 denotes the norm. The second says that under the same conditions for any 
derivative of order s 

WD^ifg) - fD'^gh < C{\\D^ fhWgWoo + ||i^/||oo||i^^-'^||2) (4.6) 

The last estimate concerns the composition of a smooth function F with a function in a 
Sobolcv space. Suppose that F is a smooth function defined on an open subset U of R'^ 
and that / takes values in an open set V of R'^ whose closure is compact and contained 
in U. Then for s > 1: 

\\D%Fif))h<C\\F\\cs\\fr^'\\D^fh (4.7) 

Here the C"* norm of F is taken over V. Note that s = is excluded in (4.7). An inequality 
for the case s — can be derived as follows. Choose some fixed uq G U. Then there exists 
a smooth matrix-valued function M onU xU such that F{u) = M{u,uo){u — uo) (cf. [11], 
p. 77). It follows in an elementary way that: 

||i^(w)||2 < \\F{uo)\\2 + \\M{u, uo)\\oo\\u - uoh 

Now, for fixed uq, ||M('u, wo) ||oo can be bounded in terms of the norm of F on F for 
in any open set V of the type introduced above that contains Hwoll- If in addition U 
contains the origin and F{0) = then the inequality reduces to ||F(t(,)||2 < ||w||2, 
which is a rather close analogue of (4.7). The fundamental tool used in the existence proof 
for the nonlinear equations is the following standard energy estimate for a linear symmetric 
hyperbolic system. 
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Lemma 4.1 Let w be a smooth solution of the hnear hyperbohc equation A^{x)dou + 
A°'{x)daU + B{x) = whose restriction u{t) to each hypersurface t =const. has compact 
support. Suppose further that the functions A°' — A" have compact support, where A"^ 
are constant matrices with A^ positive definite. Then: 

\\u{t)\\Hs < \\u{0)\\hs 

+ C [ [{\\A^\\c^ + \\dtA^\\co)\\u\\Hs + {\\Du\\oo + \\dtu\\^)\\A^\\Hs + \\B\\Hs]it')dt' 

Jo 

(4.8) 

where the constant C only depends on upper and lower bounds for the quadratic form 
defined by A°. 

Proof Applying the derivative D°' to the equation gives: 
A^doiD^'u) + A^'daiD'^u) + D'^B 

= -[D'^iA^dou) - A^D'^idou)] - - A"L>"(a„w)] (4.9) 

= say. 

This can be used to compute the time derivative of the quantity f {A^D°'u, D'^u). Differ- 
entiating under the integral and substituting in (4.9) gives a sum of terms of which only 
one contains derivatives of u of order higher than that of D". However this derivative can 
be eliminated by partial integration: 

J {-A''da{D"u),D''u) = ^ J {daA''D'^u,D''u) (4.10) 

The equation which results is: 

d/dt{ J {A^D^'u, D^'u)) = J {{dtA^ + daA'')D'^u - 2D'^B - D'^u) (4.11) 

Note that the inner product defined by A^ is uniformly equivalent to the standard inner 
product on R'^. Hence it follows that, if iV" = ^y{A^D°'u,D°'u), then 

\d/dt J (iV")2| < J [N^iWdtA^^ + a„A"||ooiV" + WB\\2 + ||Q"||2)] (4.12) 

The norm of Q°' can be estimated with the help of the Moser estimate (4.6). Using 
this fact and integrating (4.12) with repect to t gives: 

(Ar«)2(t) = / (Ar«)2(o) 

t (4-13) 

+ J {iV"(t')[||9.^''||oo||«|k« + \\d,u\\oo\\DA''\\Hs-i + ||L>"S||2]K 

Formally, d/dt{{N'^ f) = 2iV"d/dt(iV"), so that one factor A^" can be cancelled in this 
formula. This formal calculation can be justified. Adding the inequalities (4.13) for all 
derivatives D°' of order less than or equal to some fixed s completes the proof of the lemma. 
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Proof of theorem 4.1 Define uq and /o to be the time independent functions which agree 
with the initial data on the hypersurface t = 0. Then define an iteration recursively as 
follows. If Un and fn have been defined, substitute these for u and / in all places in the 
reduced Einstein equations (including the energy-momentum tensor) except those where 
derivatives of u occur. Replace these derivatives by the corresponding derivatives of ttn+i- 
This gives a linear hyperbolic equation for Un+i on an interval [0, Tn+i) which can be solved 
with the initial datum which is to be prescribed for u. Here T„+i is the largest time such 
that the metric defined by certain components of is non-degenerate. Next substitute 
Un+i into the Vlasov equation and replace / by fn+i- This gives a linear hyperbolic 
equation on the interval [0, T^+i) which can be solved with the initial datum which is to 
be prescribed for /. This defines Un+i and fn+i on the interval [0,T„_|_i). In this way 
approximate solutions (w^, fn) of the reduced system are obtained. They all induce the 
correct initial data on the initial hypersurface. Let C^([0, T], i?*(R"^)) denote the Banach 
space of functions on the interval [0, T] with values in the Sobolev space H^{IV^). It 
will be shown that if s > 5 and T > is chosen sufficiently small then Tn>T for all n, the 
sequence tt„ is bounded in the space C°([0, T], ii/'^(R^)) and the sequence fn is bounded 
in the space C*^([0, T], if^(R^)). The essential point is that for s > 5 we have the Sobolev 
embedding theorem which says that in R^ and R^ any function is and there is a 
constant C such that HwUci ^ Hence, if the estimate of Lemma 4.1 is applied to 

the equation for ttn+i, the pointwise norms of and u occurring there can be estimated 
in terms of the norms of the same quantities. To estimate the norm of diA^{un) which 
occurs, first apply the chain rule and then substitute for the term dtUn which comes up 
using the equation. Finally, apply the Sobolev embedding theorem to the result. The i?^ 
norms of the coefficients A^{un), A°-{un) and the part of B{un) which does not involve 
the matter quantities can be bounded by a polynomial in the norms of tt^, using the 
Moser estimate (4.7), on any interval where u takes values in a compact subset of U . (Note 
that this part of B maps the origin to itself and so an estimate is also obtained for the 
undifferentiated quantity, as in the discussion following (4.7).) It is also straightforward to 
show that, in the presence of a bound for the maximum momentum P{t) of any particle 
in the support of the distribution function, the norm of the energy-momentum tensor 
can be estimated by a constant depending on the norm of u times the norm of /. 
The result of all this is that if Unif) = sup^<^ ||'i^fe(^)||//^ and Fn(t) = sup;j<^ ||/fc(^)||H'^ 
then inequalities are obtained of the form 



where G and H are polynomials with non-negative coefficients. These inequalities are 
obtained under the assumption that the Uk with k < n take values in a compact subset 
K of U and that the Pk are all bounded by some constant P. The functions F and G 
may depend on K and P. The inequalities are reminiscent of (3.26) and (3.27) and once 
again the solutions of the integral inequalities can be compared with the solutions of the 
corresponding integral equations. Choose the number T that the solution of the integral 




(4.14) 
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equations with initial data exists on the interval [0,T]. For n > 1 the 

validity of these inequalities for the functions Un and depends on knowing that the Uk 
take values in a fixed compact set K ofU for all A; < n. It can be proved by induction that 
these inequalities are valid for T sufficiently small. Note first that they obviously hold for 
n = 0. Now suppose they hold up to some given value n of the index. Then they give a 
bound on ||5tt(n||oo- This in turn can be used to show that for T small enough (with a 
smallness condition which does not depend on n) Un takes values in a fixed compact set K 
and is less than a fixed constant P. This completes the inductive step. It follows 

that Un and are uniformly bounded on the interval [0,T]. 

Next note that if the coefficients A"^ of a symmetric hyperbolic system are smooth 
there exist smooth functions A" such that A'^{u) — A"{u') = Af{u, u'){u — u'). A similar 
statement applies to the function B{u). Hence the difference Wn+i ~ satisfies a linear 
symmetric hyperbolic equation whose coefficients involve A°('Un, Wn-i)- The difference 
fn+i — fn satisfies a similar equation. Using the estimate of Lemma 4.1, the uniform 
bounds for Un and Fn and the fact that the initial data are the same for all iterates leads 
to an estimate of the form: 

\\Un+l - Wn||ii-=-i + ll/n+1 " /n||i?«-i < CT{\\Un - Un-l\\H=-^ + \\fn - /n-l||j/«-i) (4.15) 

If T is chosen small enough this shows that the sequences Un and fn are Cauchy sequences 
in the Banach space C°([0, T], and hence converge to some limits {u, f) in that 

space. It should be noted that the analogue of the estimate (4.15) with H^~^ replaced 
by does not follow from Lemma 4.1. The reason is that the equation satisfied by the 
differences of iterates contain first derivatives of the iterates as inhomogeneous terms. To 
get some statements about convergence in i?*, a little functional analysis will be used. 
This uses once again the fact that the iteration is bounded in the norm. The Banach- 
Alaoglu theorem ([18], p. 115) and the fact that L°°([0, T],H'^) is the dual of a Banach space 
implies that there is a subsequence which converges weakly in L°°([0,T], H^). It can only 
converge to (w, /). Thus the limiting functions u and / are in L°°([0, T], if*). If s is chosen 
to be at least seven then Un and /„ and their spacetime derivatives converge uniformly 
and hence the limiting functions satisfy the equations. It must still be checked that if the 
initial data satisfy the condition d^gctf3 = ho-py, then the solution also has this property. 
This holds because the equations (4.3) imply the equation do{haf3'y ~ d^goip) = 0. The time 
of existence depends only on the norm and an argument used in the proof of Theorem 
2.1 shows that as long as the norm, s > 7, is bounded on some interval the solution 
can be extended to a longer time interval. If, in deriving a differential inequality for the 
solution, we keep the norms instead of eliminating them using the Sobolev embedding 
theorem, we see that the reulting inequalities are linear in the H"^ norm. What this means 
is the following. If any norm, s > 7, is bounded then the norm is bounded and then 
the norm is bounded for every s. Hence the solution can be extended in the space 
for every s. This implies that a solution corresponding to C°° initial data is itself C°° as 
long as it exists in . The statement concerning the support of / was already proved in 
Section 1 and hence the theorem is proved. 

The estimates used in this proof can also be used to prove a local existence and uniqueness 
theorem for data of finite differentiability. The idea is to approximate the data of finite 
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differentiability by a sequence of C°° data which converge to the original data in some 
Sobolev space . By Theorem 4.1 there is a C°° solution corresponding to each of these 
initial data sets and it can be assumed without loss of generality that all these solutions 
exist on a common time interval [0,T], since the sequence of data is bounded in . If 
u and u' are solutions of a symmetric hyperbolic system then their difference satisfies a 
homogeneous linear symmetric hyperbolic system. The argument is the same as that used 
to estimate the difference of iterates above. The estimate of Lemma 4.1 can be applied to 
this equation to show that if u and u' are bounded in it is possible to bound the H^~^ 
norm of u — u' by a constant multiple of the corresponding norm of the initial data. A 
similar estimate can be obtained for the reduced Einstein- Vlasov system. It follows that 
the sequence of C°° solutions introduced above converges to a solution corresponding to the 
initial data of finite differentiability. It may be remarked in passing that the uniqueness of 
the solution whose existence is asserted by Theorem 4.1 can also be proved by this method. 

It is possible to prove energy estimates in weighted Sobolev spaces with a weight 
5 > —3/2 in a way which is similar to that which has been done above for ordinary Sobolev 
spaces. The weighted Moser estimates which are necessary can be obtained using the tech- 
niques of Bartnik[2] from the Moser estimates for a bounded region in R"^. These estimates 
imply the boundedness of a sequence of approximating solutions in L°°([0, T], Hf). Apply- 
ing the Banach-Alaoglu theorem once more shows that the previously obtained solution is 
in this space if the data are in H^. This is still not the most precise result on the prop- 
agation of aymptotic fiatness which is desirable. For it is usually assumed that in initial 
data sets for the Einstein equations, the second fundamental form falls off faster than the 
spatial metric and it is desirable that this property be preserved by the time evolution. 
For a general symmetric hyperbolic system, there is no obvious reason why this should 
be true. For if one tries to use the equation directly to prove a property of this kind, 
the term B{u) intervenes and in general there is no reason that this should fall off faster 
than u itself. However if B at least quadratic it is true. Here this is only sketched briefly 
and we do not even give a precise definition of the phrase 'at least quadratic'. If B is at 
least quadratic and u G and 5 > —3/2, then B{u) belongs to Hg, for some 6' > S. 
If 5' < 5 + 1 then it follows from the equation that dtu is in Hg~^. If S' > 5 + 1 then 
dtu is in Hg~l. If necessary this procedure can be repeated and after finitely many steps 
the result is obtained that dtU G Hg~l. In the case of the Einstein- Vlasov system, the 
matter term in the Einstein equations is irrelevant for this discussion since its support is 
in a known compact set. On the other hand, the part of B only involving the geometry is 
quadratic in the Christoffel symbols. The energy estimates show that the time of existence 
of the solution in a weighted Sobolev space only depends on a weighted Sobolev norm of 
the initial data. 



Appendix 

Lemma Al Let / : [0, oo) ^ R be a function, q a polynomial on R"^ which does not 
vanish identically and is homogeneous of degree 2m, where m is an integer, and let F be 
the function on R" defined by F{x) = f{\x\)q{x)/\x\^'^. If m > 1 suppose further that / 
is (7^*^"^ and /(^'^(O) = for aU integers I with I < m. Then for any A; > the following 
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conditions are equivalent: 

(i) / is and = for all integers I such that 21 + 1 < k 

(ii) F is 

Proof If (ii) holds consider g{x^) = F{x\0,...,0) = /(|x^|)[g(x\ 0, . . . , 0)/(x^)2"^]. The 
second factor is a constant. If this constant is non-zero it can be seen immediately that / 
is C^. If the constant is zero it can be made non-zero by a rotation. Since also g{x^) = 
g{—x^), the statement (i) holds. Conversely, suppose that (i) holds. Let pk be the Taylor 
polynomial of order A; of / at the origin. Then we can write: 

Fix) = pki\xMx)/\x\^-^ + ifi\x\) - pki\x\Mx)/\x\'-^ (Al) 

Taylor's theorem applied to / shows that the second term is with all its derivatives 
up to order k vanishing at the origin. To complete the proof it only remains to show that 
under the hypothesis (i) the first term is a polynomial. Under this hypothesis the first 
2m — 1 derivatives of pk at zero vanish, as well as higher derivatives of odd order. Hence 
Pfc(|a;|) is a polynomial in |xp which is divisible by and the first term in Al is a 

polynomial. 

Lemma A2 Suppose that the hypotheses of Lemma Al are satisfied and that the deriva- 
tives of / of order < k are bounded. Then the partial derivatives of F of order < k are 
bounded and there exists a constant C independent of / such that ||-F||cfc ^ ll/llc*^- 
Proof On the region |a;| > 1 the statement is obvious. To see that it holds when r < 1 
is included, first split F into a sum of two terms as in Al above. The first term and 
its derivatives can easily be bounded in terms of derivatives of /. A derivative of the 
second term of order < /c is a sum of terms, each of which is the product of a derivative 
of if — Pk){\x\) with a function which is homogeneous of degree > —k and is independent 
of /. The derivative of (/ — pk){\x\) of order I < k can be bounded on the region |a;| < 1 
in terms of the |a;|'. The resulting powers of \x\ are sufficient to compensate the 

functions which are homogeneous of negative degree and this gives the desired estimate. 

Lemma A3 Let / : [0, oo) — > R be a function with /(O) = and let be the vector field 
on defined by F'^ix) = fdxDx" /\x\ for a; and F(0) = 0. Then for any /c > the 
following conditions are equivalent: 

(i) / is and f'^^\0) = for aU I such that 21 < k 

(ii) F« is C'' 

Proof F^{x^ ,0, . . . ,0) = f{x^) for x^ > and this immediately shows that (ii) implies 
that / is C^. Then the relation F^{—x^) = —F^{x^) shows that the derivatives of even 
order vanish, proving (i). Conversely, suppose that (i) holds. The vector field F"' can be 
written as a sum of two terms, as in the proof of Lemma Al, and the second of these terms 
can be handled just as in the proof of that lemma. As for the first term, it is the product 
of a;°'/|a;| with a polynomial in |a;| where only the coefficients of odd powers are non-zero. 
Hence it can be written in the form x°'q{\x\^), for some polynomial q. 

Lemma A4 Suppose that the hypotheses of Lemma A3 are satisfied and that the deriva- 
tives of / of order < k are bounded. Then the partial derivatives of F'* of order < k are 
bounded and there exists a constant C independent of / such that ||F"||crfe < H/Hc*:- 
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Proof The method used in the proof of Lemma A2 also apphes in this case. 
Lemma A5 Let / : [0, oo) — > R be a measurable function and define 

Jo 

If w > —1, V < w and / is then g^^^ is C*^"*"-*^. Moreover, on any interval of the form 
[0, R] there exists a constant C, independent of /, such that q||(^fc+i < C||/||cfc. 
Proof Write f — Pk + if — Pk), where pk is the Taylor polynomial of / of order k at the 
origin. Now r~'" s^pk{s)ds can be bounded by a constant times ll/Hc* under the given 
assumptions. Thus it can be assumed that without loss of generality that pk — 0. Now: 

d^9v,w/dx^ = d\r-'')/dx^ [ s"" f{s)ds + y^hjd^-^ f /dx^'^ 

Jo 

where hj is homogeneous of degree —v + w + j. Under the assumption that pk = the 
integral in the first term is bounded by Cr^'^'^'^^, where the constant C only depends 
on the norm of /. Since the other factor is 0{r~'"~^) this suffices to bound the first 
term. The other terms can be bounded using the fact that d^~^ f /dx^~^ is bounded by an 
expression of the form C\x\'^~^'^'^. 

Lemma A6 Let / : [0, oo) — > R be a continuous function with /(O) = and define 

h{r) =r [ s-^f{s)ds 
Jo 

If / is for some k > 1 then h is C^^'^. Moreover, on any interval of the form [0,R] 
there exists a constant C, independent of /, such that ||/i||(^fe+i < C||/||cfe- 
Proof Since /(O) = and / is for some /c > 1 it is possible to write /(r) = r/i(r) for 
some C'^~^ function fi. Then h{r) = r fi{s)ds. Clearly r~^h is C'' and its C'' norm 
can be estimated by the C'^ norm of /. It follows that h itself is C'=+^ and that its C'=+^ 
norm can be bounded by the norm of /. 
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